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1. Introduction 

T-H ■ Recent observations of the corona with the LASCO coronagraph 

[ (Schwenn et al., 1997) on board of the SOHO spacecraft showed that 

^ ■ the corona can be highly structured even during the solar activity 

, minimum. The observations revealed a triple structure of the streamer 

CN I belt which was existent for several consecutive days. The observations 

^ ■ further showed that these triple structures occasionally go unstable 

Q>^ i leading to a seemingly new and extraordinarily huge kind of coronal 

' mass ejection (global CMEs). Natural questions arising from these ob- 

Q^. servations are whether the helmet streamer triple structure is directly 

^ I connected with or responsible for the occurrence of global CMEs and 

^ ■ what is the physical mechanism of their formation. The aim of this pa- 

"to . per is to provide a first step towards a better theoretical understanding 

' of these phenomena. 

^ ■ The structure of helmet streamers and their stability has been 

^ I studied both observationally and theoretically for a long time (e.g. 

^ • Pneuman and Kopp, 1971; Cuperman et al., 1990; Cuperman et al., 

1992; Koutchmy and Livshits, 1992; Wang et al., 1993; Cuperman et 
al., 1995; Wu et al., 1995; Hundhausen, 1995; Bavassano et al., 1997). 
A natural association seems to exist between helmet streamer stability 
and coronal mass ejections (CMEs) and sometimes CMEs are accompa- 
nied by erupting prominences or filaments. Coronal streamers are also 
thought to be the source regions of the slow solar wind and the activity 
processes discussed above may well contribute to this component of the 
solar wind. 

Multiple streamer structures have already been observed before 
SOHO, however mainly during the maximum phase of the solar activ- 
ity cycle. Also, observations of the heliospheric plasma sheet seem to 
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indicate a multiple current sheet substructure of the plasma sheet itself 
(Crooker et al., 1993; Woo et al., 1995) which has initiated several the- 
oretical investigations of the stability of multiple current sheets (Otto 
and Birk, 1992; Yan et al., 1994; Dahlburg and Karpen, 1995; Birk and 
Otto, 1997; Birk et al., 1997; Wang et al., 1997). 

The aim of this work is to present a simple method by which triple 
helmet streamer structures can be calculated self-consistently. We want 
to keep the method simple enough so that we can obtain largely ana- 
lytical results without rendering our models too unrealistic. Since the 
observed structures seem to be quite elongated (the radial magnetic 
field component is considerably larger than the longitudinal magnetic 
field component), we can employ a method which has already been 
applied very successfully to the equilibrium structure of the magnetotail 
of the Earth (Schindler, 1972; Birn et al., 1975). The method is based 
on an asymptotic expansion. Here we generalize this method and show 
how it can be applied to find solutions with internal multiple structures. 
Obvious generalizations of the original magnetospheric versions to solar 
applications are the use of a more suitable coordinate system (spherical 
rather than Cartesian), the inclusion of magnetic shear, of the solar 
gravity field for large scale structures and of plasma flow on open field 
lines. For reasons of mathematical simplicity and since we only deal 
with the subcritical part of the flow, we confine the discussion to incom- 
pressible flow. Our models provide a flexible quantitative description 
of structured helmet streamers. They can also be used as the starting 
configurations for future studies of the dynamical evolution of such 
structures. 

The outline of the paper is as follows. In Sect. ^ we present the math- 
ematical basis of the method. Section |^ contains representative results 
obtained by applying the method under various sets of assumptions. In 
Sect. Q we summarize the paper, give a discussion of our results and an 
outlook on future work. 



2. Mathematical Formulation 
2.1. Basic equations 

We use the equations of stationary ideal magnetohydro dynamics to 
describe the coronal plasma: 

- VP + — (V X B) X B - pV^' = pv • Vv (1) 

Mo 

V • (pv) = (2) 

E + V X B = (3) 
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V-B = 



P = pRT. 



(4) 
(5) 



Here, P stands for the plasma pressure, B for the magnetic field, p 
for the plasma density, ^ for the solar gravitational potential, v for 
the plasma velocity, E for the electric field, R is the gas constant, 
T the temperature and po the vacuum permeability. Due to the high 
conductivity of the coronal plasma the magnetic field is frozen into it 
(Eq. ^) and we can make the assumption that the observed plasma 
features outline the magnetic field structure. Throughout the rest of 
the paper we make the following assumptions: 

— The streamer structure is very extended in azimuth and varies 
only very slowly along this direction (this assumption is supported 
by the observations (Schwenn et al., 1997) which showed that the 
triple structure existed for several days); we can then describe 
the structure as approximately two-dimensional in the sense of 
rotational invariance. 

— The streamers are elongated in the radial direction and the radial 
magnetic field component is considerably stronger than the lati- 
tudinal component; this is the basic assumption allowing to apply 
an asymptotic expansion procedure. 

— If flow is included, it is purely field-aligned and confined to the 
open field line region; since we describe only the low lying parts 
of the streamer structure the flow velocities should be subcritical 
and we will therefore only consider incompressible flow. 

— For simplicity, we consider only isothermal equilibria; in principle, 
other temperature structures or even the inclusion of an energy 
equation into the scheme is possible. 

In the following, we normalize the magnetic field by a typical value Bq, 
the plasma pressure p by Bq / pQ (where p = 1 is the maximum value 
of p), the mass density p by Bq/pqRT, the length L by a solar radius 
and the current density by jy = Bq/pq L. 

2.2. Method 

To demonstrate the method we first use a two-dimensional Cartesian 
geometry with no spatial variation in the y-direction (^ = 0) and 
assume static structures (v = 0) without magnetic shear {By = 0). In 
other geometries the basic properties of the method stay similar though 
the details may vary due to the geometry. 
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Writing the magnetic field as 

B{x,z) = VA{x,z) X ey (6) 

we find from equation (||) that the flux function A has to obey the 
equation (e.g. Low 1975; Birn et al. 1978; Low 1980; Birn and Schindler 
1981) 

-AA=-^iP{A,^)) (7) 

With the assumption of a constant temperature and because = '^(z) 
in Cartesian geometry we get: 

PiA,"^) = P{A,z) = k{z)p{A) (8) 

with k{z) = exp(— (^'(z) — ^'(0))/RT). For slow variation of ^{z) with 
height z one can neglect gravity and finds k{z) = 1. 

If we assume that the magnetic field perpendicular to the pho- 
tosphere {Bz in the Cartesian geometry) is much larger than the 
parallel component (B^ in the Cartesian geometry), or in other 
words, that streamers are rather elongated in the radial direction, 
we can use the method of asymptotic expansion and apply it to Eq. 
(^). Mathematically stretched configurations are characterized by the 
ordering 

_ = 0(1) ^B, = 0(1) — = 0(e) « 1 ^ i?, = 0(e) (9) 

If we neglect terms of the order e^ in Eq. we obtain after one 
integration with respect to x : 

^ = ±y^2{po-p{A)) (10) 

where pq = po{z) is an integration constant (total pressure on the z- 
axis), which depends parametrically on z. We solve this differential 
equation by separation (Birn et al., 1975): 

dA 

X — Xq = / — , = (11) 

Jao y'2{po-p{A)) 

and fix the constant = ^(^^o) so that p{Aq) = pQ. 

Unfortunately it is impossible to derive the pressure function p{A) 
directly from the observations presently available. We will therefore 
choose p{A) by a compromise between physical reasoning and mathe- 
matical simplicity. We do not choose pressure functions which have a 
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too simple dependence on A like p{A) oc const., A or A"^ because solu- 
tions of this type are a) known to be linearly stable and b) have already 
been discussed in Wiegelmann (1997). A popular and convenient choice 
is p{A) oc exp(c^) because it corresponds to the physical assumption 
that the plasma is in local thermodynamic equilibrium and it has the 
mathematical advantage that one can find analytical solutions. Using 



p{A) = k exp(cA) 



fl2) 



where k may be a function of z and c may be positive or negative, we 
get from equation (|lT|): 



A{x, z) = — log < cosh 



+ - log ^ 



c\/ - ^o) 



The components of the magnetic field are then given by 



1 



Po 



dA 



dx 



dA dA dpo dA dk 

dz 



(13) 



(14) 



(15) 



dpo dz dk dz 

We remark that one cannot evaluate and ^ in equation (15) 
directly because in general xq depends on po{z) and k{z) as well. The 
plasma pressure can be obtained from the Eq. ( [T^ ) 



Po 



P = , 

cosh2( c(x - xq)) 

and the current density jy from jy = ^j- 

Poc 



(16) 



3y = cp 



cosh 



-C{X - Xq) 



(17) 



2.3. Triple structures 

The method just outlined has been applied very successfully to cal- 
culate the structure of the Earth's magnetotail (Birn et al., 1975; 
Schindler und Birn, 1982; Wiegelmann and Schindler, 1995). However, 
since we want to apply this method to describe multiple streamer struc- 
tures a few details of the method have to be changed though the basic 
approach stays the same. 
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In an arcade structure the direction of the component of the 

magnetic field changes its sign as one crosses the centre of the arcade. In 
a triple structure, one will therefore encounter three such changes as one 
moves across the structure. In a stretched configuration described by 
the asymptotic expansion procedure, the only possibility to achieve this 
is that the current density jy changes its sign between the streamers, or, 
in other words, in the central streamer the current flows in the opposite 
direction with respect to the two outer streamers. In terms of the sign 
of Bz this means that in order to paste the three arcades continuously 
together, the B^ component in the central streamer has to change just 
in the opposite sense to the outer two streamers. 

Since the current density jy is directly linked to the plasma pressure 
by the relation jy = dP/dA, this implies that the derivative of the 
pressure with respect to the flux function A has to change its sign, too. 
We can therefore not use a single exponential function to model the 
pressure because the exponential function and its derivative have no 
zero. The simplest alternative is that we use different pressure func- 
tions for the middle {pi{A) say) and for the outer streamers {p2{A)). 
In principle, we could also use a different pressure function for each 
of the three streamers, but since we do not want to make the model 
too complicated we assume the structure to be mirror symmetric with 
respect to the z-axis for simplicity. 

Due to the different pressure functions our solution will have two 
separatrix fleld lines characterized hy A = Ag which separate the re- 
gions of different pressure from each other. These separatrices have to 
be calculated self-consistently together with the solution but due to our 
symmetry assumption we only have to calculate one of those separatrix 
field lines. In the asymptotic expansion procedure the location of the 
separatrix field line will be given by Xgep = /(z) where / is only varying 
weakly with z. 

If we require a smooth transition of magnetic field across the separa- 
trix from the outer to the inner streamer, the pressure functions must 
be continuous. Therefore we choose 

pi(^) = A:iexp(-ci^) (18) 

for the middle streamer and 

P2{A) = k2cxp{c2A) (19) 

for the outer streamers, while the coefficients ki, k2, c\ and C2 (ci, C2 > 
0) are related by 

k2 = kx- exp(-(ci -I- C2)As) (20) 

on the separatrix labeled A = Ag. This choice has the consequence 
that the current density jy has a discontinuity at the separatrix. We 
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Figure 1. This picture illustrates the right side of our configuration (the left hand 
side can be derived from the symmetry assumption). The dashed ( ) line cor- 
responds to the separatrix Xsep{z) between the streamers, the dotted (■ ■ ■) line 
corresponds to the center of the outer streamer 3:02(2) and the dashdotted (— •) 
line corresponds to the outer separatrix Xout{z). 



remark, however, that the absolute value of jy is only rather small at the 
separatrix and only changes from a small negative to a small positive 
value and therefore the discontinuity of the current density should not 
influence the solution structure very much. 

We choose xq = to put the centre of the middle streamer on the 
z-axis. To calculate the separatrix between the middle and one of the 
outer streamer we use equation ( [TTl) with the form of p{A) inside the 
middle streamer to obtain 

larctanh ^^^Ih^t^^ . (21) 

Note that both po and ki may depend weakly on z. 

If we call the centre of the outer streamer X02, we can calculate this 
central field line relative to the separatrix field line again with equation 
(|lTl), but this time using P2iA) and setting xq = xo2- 

X02 = Xsep + \ ——Si^'cctanli i \ — — ^2 exp(c2^s) | ^22) 
V Po C2 y\j PO J 
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If one takes into consideration that pi{As) = p2{As) and Equation (21) 
one finds: 



= arctanh (P. - fa exp(-c.^) /I 1 + i , (23) 



arctanh 



{po - k2exp{c2As) 

PO 




A sketch of the whole structure is shown in Figure 0. 

We can now calculate A, B^, Bx, p and jy as functions of x and 
z from the Equations (jl|), (0), (|l5|), (|l|) and (|l^ and by using the 
symmetry properties of these quantities. As an input we need the total 
pressure po{z) which is the same in all three streamers, the parameters 
ci, C2 and ki and k2- If gravitation is included, ki and /c2 may depend 
weakly on z, but could be different if we allow different temperatures 
inside the streamers. 



2.4. Inclusion of flow 



The method as outlined so far neglects that coronal streamers are 
the source regions of the slow solar wind. It is an open question if 
a stationary slow solar wind component exists or whether the slow 
solar wind is the result of many small eruptions. Recent observations 
(Schwenn and Inhester, private communication 1997) seem to favour 
the second alternative. If that is true, a steady state model providing 
the start configuration for a time-dependent computation of the erup- 
tion process need not include the flow, at least not near the symmetry 
plane. Nevertheless we want to show how a stationary plasma flow can 
be included in our model. It may be relevant for the open field regions 
at larger distances from the ecliptic, to which we confine the flow. 

The inclusion of plasma flow into the asymptotic expansion proce- 
dure has been carried out for both incompressible (Gebhardt, 1989; 
Birn, 1991) and compressible flow (Birn, 1992; Young and Hameiri, 
1992). There are three critical points in the solar wind: The slow 
magnetosonic point Vg, the Alfven point va and the fast magnetosonic 
point Vf with Vs < va < Vf for the coronal plasma. If one assumes 
that the flow velocity at the basis of the corona is less than Vs the solar 
wind has to pass through all these critical points. Birn (1992) showed 
that in stretched configurations the total pressure function po{z) is 
related to the flow velocity. For v < Vg the flow velocity decreases 
with decreasing pq{z) and for v > Vg the flow velocity increases with 
decreasing Pq{z). As it is necessary that po{z) decreases with respect 
to z to get closed arcade structures, the solar wind velocity should 
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have a minimum at the helmet streamer cusp. This result is of course 
only valid if a stationary slow solar wind exists. The assumption of 
incompressible flow in the solar wind, however, cannot be justified 
in a strict sense, because generally gravity is not negligible. Since in 
the paper our main emphasis is placed on closed structures without 
directed flow, this model assumption does not influence our conclusions 
significantly. (In fact, all our explicit examples are structures without 
flow.) We note that the above mentioned variation of the solar wind 
velocity remains qualitative valid, if one neglects gravity and assumes 
incompressible flow. We will assume that the flow is purely field-aligned 
and since we will not extend our solution to the first critical point, we 
will for simplicity make the assumption of the flow being incompressible 
and subalfvenic. 

We mention that, as a useful nice property, such stationary solutions 
with flow, can be directly calculated by transformation from static 
solutions (Gebhardt and Kiessling, 1992). Here we follow the approach 
developed in Birn (1991). Equation (|lO|) is then substituted by: 



dA_ l 2{po-p{A)) 

dx~ \l 1 - MliA) ^ > 

where Ma{A) is the Alfven Mach number defined by M4 = v/va, where 
VA = B/ ^/Jiop is the Alfven velocity. Note that for incompressible flow 
both the Alfven Mach number and the plasma density are constant 
along open field lines. 

To extend the method outlined so far to include flow on open field 
lines, we have to define a second separatrix field line between the 
regions of open and closed magnetic flux. Since there shall be no flow in 
the closed field line regions, we choose Ma{A) = on the left hand side 
of the separatrix and < Ma{A) < 1 (subalfvenic flow) on the right 
hand of the separatrix. In the special case Ma = constant we find with 
equation (^) that one can easily include incompressible subalfvenic 
plasma flow (Ma < 1), if we replace the term {x — xq) in the static 



theory by {x — xq)/^ 1 — M^. The location of the separatrix field line 
Ay could be calculated by a similar procedure as the inner separatrix, 
but it is more convenient to choose it to be at the same distance from 
the centre of the outer streamer as the inner separatrix so that we have 
a symmetric situation again and one finds A^ = As- We remark that 
this condition is necessary if one investigates configurations with a cusp 
structure. One can then calculate the shape of the separatrix defined 
by ^outiz) simply from XQ^ti^z ) = 2xq2{z) — Xsep{z) (see Figure |l] for 
illustration) . 
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2.5. Inclusion of magnetic shear 

It is possible to extend our method and include magnetic shear writing 
B as: 

'B{x,z) =VA{x,z) X Cy + By{x,z)ey (25) 

Using the y component of equation (||) one finds By = By{A{x, z)) and 
instead of (^) we get: 

-AA = -^iU{A,^)) (26) 
with ^ 

n(A, ^) = U{A, z) = k{z)p{A) + ^4^-- 
Using the method of asymptotic expansion we get: 

_ dA 

Jao y/2{po{z)-U{A,z)) 

We remark that the reason for magnetic shear should be a displacement 
of magnetic foot points and that only closed field lines will be sheared. 
We assume By{A) oc exp(d^) for mathematical simplicity. In general 
we cannot solve equation (|2^) analytically. This is only possible if we 
choose the special case d = § which we will do to illustrate the effect 
of magnetic shear. Then we get: 

n(A^) = (fc(^) + y)exp(cA) (28) 

and we can easily include magnetic shear by substituting k{z) by k{z) + 



^ in equation ([131) . It is possible to investigate other values (c ^ 2d) 
numerically, but one can see the main effect of magnetic shear already 
in this special case. 

We calculate triple structures analogous to the method without 
shear and substitute p{A,z) by 11(^,2:). We have to consider shear 

only on closed field lines. We use ni(yl) = (^k{z) + exp(— ci^) in 

the middle streamer, n2(^) = exp(— (ci +C2)As) (j^iz) + exp(c2j4) 

in the outer streamers and n3(A) = exp(— (ci + C2)As)/c(z) exp(c2A) on 
open field lines. Thus there is a jump in 11 and a corresponding jump in 
the magnetic field at the boundary between outer streamers and open 
field lines and consequently there is a thin current sheet. Unlike the 
calculations on closed field lines one cannot fix U^^Aq) = pQ. Instead 



of (11) one has to use: 



xo= [ ^ (29) 

J A. V2(po-n3(^)) 
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where is the value of the last closed field line (boundary field line 
between open and closed regions). One finds = Ag for configurations 
with cusp structure. 

We remark that violating the condition ni(j4s) = 112 (^4^) could 
still make sense. For example, one could use Pi{As) = P2{As) but 
Byi{As) 7^ By2{As). Of course this mismatch in the value of IT would 
have to be compensated by a corresponding mismatch in Bz and this 
would lead to additional current sheets between the middle and the 
outer streamers. 



2.6. Spherical coordinates 



In this subsection we formulate the method in spherical coordinates 
(r, which are more realistic to describe the solar corona on scales 
comparable to or larger than a solar radius. We assume static (v = 0) 
configurations and do not consider spatial variation in the (^-direction 
= 0). Thus we can present the magnetic field as: 



Br 



B, 



1 

sin 6 86 
1 



r sin 6 dr 

B^ = B4r,d) (30) 

Here ^ is the flux function in spherical coordinates. We remark that in 
spherical coordinates the flux function is not identical with a compo- 
nent of the magnetic vector potential. One finds = r sm(0)A^, where 
A^ is the (/>-component of the vector potential A. We find from ([l|) that 
^(r, 6*) has to obey the equation: 

dr^ r2 86 \sme 86 J 3^ 2 8^ ^ ^ 

where /(^) = rsin^i?^ is constant on field lines. With the assumption 
of a constant temperature and because ^ = ^(r) one finds: 

P{C,^) = k{r)p{0 (32) 

with k{r) = exp(-(^'(r) - *(0))/RT). 

We assume that the magnetic field perpendicular to the photosphere 
Br is much larger than the parallel component Bg. For simplicity we 
do not consider magnetic shear here {Bt^ = 0) . Stretched configurations 
are characterized by 



\§0 = ^Br = 0(1) i|- = 0(e) « 1 ^ 5, = 0(6) 
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This approximation is valid for configurations close to the equatorial 
plane {6 = ^) and r'^po{r) has to vary only slowly with respect to r. The 
second condition limits the outer boundary of our model corona to a 
few solar radii. As the observed configurations are indeed concentrated 
near the equatorial plane and the closed streamer field lines seem to 
be radially bounded to a few (about 2 — 4) solar radii, our asymptotic 
model should describe these configurations approximately correct. 

If we neglect terms of the order in Eq. (^) we obtain after one 
integration with respect to 9: 

II =1 r2 sm{9) V2{po{r) - k{r)p{0) (33) 

We use (analogous to our calculations in Cartesian geometry) p(^) = 
exp(c^) and solve this differential equation by separation. Instead of 
LSl) we get: 



log{p(){r)) — log(fc(r) — 2 log(cosh(cA /^r^ (cos(0) — cos(0o))) 
e(r, 9) = ^ 

(34) 

We remark that one cannot find a one dimensional solution of equa- 
tion ( |3T|) with ^ = C{9) like in Cartesian geometry. Thus no analogy to 
the Harris sheet in Cartesian geometry exists in spherical coordinates. 

To calculate triple structures in spherical geometry one uses a similar 
procedure as described for Cartesian geometry. We choose 9q = ^ ^ 
sin^o = to put the centre of the middle streamer and calculate the 
separatrix field line analogous to ( ^T[) and get: 



cos(e,ep) = r%/llarctanh f /(Po^^Wz^ ) (35) 

V Po ci \^ y Po J 

If we call the center of the outer streamer 6q2 we get analogous to our 
calculations in Cartesian geometry: 



cos(0O2) = T^'^xl — ( 1 ) arctanh 

Po \Ci C2 



i + 1 I arctanh 




'{PO 


— ki exp(- 


-ciAs) 


Po 


kpo 


- /c2 exp(- 


-C2As) 



Po \Cl C2j \ V Po 



The quantities ^, Bj.-, Bg, p, can then be calculated as functions of r 
and 9 in the usual way. 
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Table I. Parameter sets used to calculate the 
solutions presented in Figure H 



Solution 


Sl 


S2 


S3 


Cl 


C2 


A, 


a) 


0.8 


0.4 


0.2 


15.0 


15.0 


0.1 


b) 


0.8 


0.4 


0.2 


15.0 


15.0 


0.08 


c) 


0.8 


0.4 


0.2 


20.0 


20.0 


0.12 


d) 


0.8 


0.4 


0.2 


10.0 


30.0 


0.12 


e) 


0.8 


0.4 


0.2 


30.0 


12.0 


0.04 


f) 


0.8 


1.2 


0.2 


15.0 


15.0 


0.12 



3. Results 

3.1. Static solutions in Cartesian geometry without 
gravitation 

We first give a few examples for triple streamers calculated with our 
model to illustrate the influence of the different functions and param- 
eters on the structure of the solutions. To calculate a solution we have 
to specify the functions Pq{z) and ki{z) and the parameters ci, C2 and 
Ag. In the first examples we neglect the effect of plasma flow and of 
the solar gravitation which gives ki{z) = 1. Our solutions are confined 
within the boundaries x = — 0.5..0.5 and z = 0..5 corresponding to one 
solar radius in the x direction and five solar radii in the z direction. 
We use symmetry with respect to the z axis. 

We prescribe the total pressure by po[z) = si exp(— S22:) + S3. Within 
our model A{x,z) depends only parametrically on z by Pq{z) and so 
the exact form oIpq{z) has not too much influence on the configuration 
as long as its general properties are similar. One important property 
of Pq{z) is that it has to be a monotonically decreasing function of z 
to obtain singly connected closed field line regions. Physically one can 
identify the total pressure on the x-axis with the sum of the magnetic 
pressure on open field lines outside the helmet streamer and a 

homogeneous plasma pressure of the solar wind. 

In Figure |2| we show field line plots of six solutions calculated for 
different sets of parameters which are summarized in Table |. In Figure 
^) - |2|c) we investigate the influence of the location of the separatrix 
field line on the solutions structure, whereas pq is kept fixed and the 
parameters ci and C2 are fixed to a value of 15.0 in a) and b) whereas 
we used ci = C2 = 20.0 in c) which leads to thinner configurations. In 
these cases all three streamers have exactly the same width though the 
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Table II. Parameter sets used to calculate the solutions 



presented 


in Fig 


;ure |3| 












Solution 


a 


b 


n 


•2^cusp 


Cl 


C2 


As 


a) 


0.8 


0.2 


3.0 


4.0 


15.0 


15.0 


0.1073 


b) 


0.8 


0.2 


4.0 


4.0 


10.0 


25.0 


0.1609 


c) 


0.8 


0.2 


6.0 


4.0 


25.0 


12.0 


0.0644 



width itself varies. In Figure |2|a) the separatrix is a straight Une giving 
rise to parallel streamers. In Figure |2|b) the separatrix has been moved 
closer to the centre of the structure and the three streamers converge, 
whereas in Figure |2|c) the separatrix has been moved further out giving 
rise to diverging streamers. 

In Figure §d) we have decreased the value of ci to 10.0 and increased 
the value of C2 to 30.0. Also the value for Ag has been increased. The 
effect is that the central streamer now contains more magnetic flux 
than the outer streamers and is wider in the x-direction. In Figure ^) 
this effect has been reversed by making C2 smaller than ci . It is obvious 
that now the outer streamers are wider than the central streamer. 

Finally, in Figure |2|f) we have varied the total pressure po{z) by 
increasing the parameter S2 to 1.2 resulting in a faster decrease of pq 
with z. One notices that the closed field line regions (which we define 
here loosely as those field lines not crossing the upper boundary) has 
shrunk and that the whole structure has a much slimmer appearance 
for large z. 

In a next step we use our model to describe helmet streamer config- 
urations with a cusp. Up to now we have not given a radial boundary 
of the closed magnetic field lines. Now we define such a point at x = 0, 
z = -Zcusp which defines the last closed field line. 

We use as a function for the total pressure: 

+ b for z < z^^^ 



6 for z > Zcusp 

Field line plots for three different sets of parameters are shown in Figure 
^. The parameter sets used are listed in Table ||. It can be seen that 
the effect of changing the parameters ci and C2 stays the same as in the 
case without cusp. An increase of the exponent n leads to a stronger 
decrease of po{z) as the cusp point is approached from below and this 
leads to a sharper cusp structure. 
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a) b) c) 




-0.5 0.5 -0.5 0.5 -0.5 0.5 



Figure 3. Magnetic field fines (contourplot of A{x, z)) for triple helmet streamers 
with cusp (see text for discussion). 

These examples illustrate the flexibility of our model. One can de- 
scribe different types of triple helmet streamers by specifying the free 
parameters in the general solutions. 



3.2. Static solutions in Cartesian Geometry Including 
Magnetic Shear and Gravitation 

In this section we investigate the influence of the solar gravitational 
field on the structure of the solutions. Though wc arc still calculating 
in Cartesian geometry we will not take the gravitational force to be 
constant as it is usually done in Cartesian geometry if the length 
scale under consideration is much smaller than a solar radius (e.g. 
Zwingmann 1987; Piatt and Neukirch 1994), but use the more appro- 
priate potential ^ = ^{z) = — f^^ ■ Together with the assumption of 
constant temperature we get: 

M.)=exp(-i-).exp(-J^) (38) 

where G is the constant of gravitation, Mg the solar mass, Rg the solar 
radius and I = . Under the assumption of a pure hydrogen plasma 
and a constant temperature of 3 • 10^ K one finds / ~ 7.5 corresponding 
to l.hRs- One difficulty of the method in Cartesian geometry is that 
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a) 1=0 1 = 



b) 1=0 X = 



c) 1=0 1=10 




Figure 4- Magnetic field fines (contourplot of A{x,z)) for triple helmet streamers. 
In the upper pictures b) and c) we included magnetic shear and in the lower pictures 
we additionally investigated the effect of gravitation. 
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k{z) must decrease more slowly with respect to z than po{z) because 
otherwise one cannot find arcade type solutions. This is difficult to 
achieve with a realistic value of I, but since we only want to investigate 
the qualitative effect of gravitation on the solution we use the unrealis- 
tic but mathematically more convenient value I = 2.5. We remark that 
this problem is alleviated in spherical geometry which is more realistic 



anyway (see section p. 31 ) 

The inclusion of shear only makes sense on closed field lines. As we 
have a well defined boundary between open and closed regions only in 
solutions with cusp structure, we only investigate such configurations 
here. In Figure ^ we fix the parameters ci = C2 
0.8, n = 4.0, b = 0.2 and use for the total pressure: 



25.0, Zcusp = 4.0, a 



k{z) + f 



-^cusp Z 
■2cusp 

h{k{z)+^ 



for 
for 



Z < Zr 



Z > Z,: 



(39) 



The factor |^A;(z) + ^ j can be used to fix the normal component of the 

magnetic field at the solar surface, which we could use as a boundary 
condition. 

The main effect of magnetic shear is that there is a jump in the 
magnetic field and consequently a thin current sheet between the closed 
and open field lines. Indications of such a current sheet are also found 
in recent observations (Schwenn et al., 1997). In Figure ^a) we plot 
as a function of x at the height z = 1 for the same parameters as in 
Figures ^) and c) . 
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If we do not include shear the structure above the helmet streamer 
cusp z > Zcusp becomes equivalent to a Harris sheet and goes 
smoothly through zero in the center of the configuration. The inclusion 
of shear on closed field lines leads to a jump in from a positive to a 
negative value in the center of the whole configuration (see figure |5|b)). 
This causes a thin current sheet, which is more pronounced than the 
current sheet at the outer boundary of the streamer configuration below 
the cusp. This current sheet corresponds to the heliospheric current 
sheet. The effect can clearly be seen in Figure § In the Figures |^) 
and c) we included magnetic shear and in c) the shear is stronger 
than in b) leading to a larger jump in Bz- This can be seen in the 
contourplot of A as the field line density is much higher in the center 
above the cusp in Figure ^) than in b). In the lower panel (Figure §d)- 
f)) we investigated the same configurations as in the upper pictures 
but included gravitation. One can see that without shear (Figure ^)), 
the configuration gets wider with increasing z. (If one uses the realistic 
value I = 7.5 instead of I = 2.5 the configuration will become unreal- 
istically wide without shear.) The inclusion of shear reduces this effect 

and if one investigates force free configurations ^^^^ = VP = 0^ or 

^) wl 

c) and f) the effect of gravitation vanishes. 
3.3. Solutions in spherical geometry 

Now we present some example solutions of helmet streamer configura- 
tions in spherical coordinates. This geometry is of course more realistic 
to describe the coronal magnetic field on large scales. In this section we 
also include the influence of the solar gravitational field on the structure 
of the solutions from the beginning. With ^ = ^(r) = —^Mm. ggf; 



nearly force free solutions {^—^^^ ^ ^^SA^ ^ which is accomplished in 



as in Cartesian geometry (see Section 3.2 



k{r) = exp y-^j ■ exp i^-j (40) 

Here we use the realistic value / = 7.5. 

In Figure ^) we show a field line plot of a solution in which we 
prescribed the total pressure as po{r) = k{r) [si ex.p{—S2{r — 1)) + S3] 
and used the same parameters as in Figure ^a) (see Table |, case a). 
Thus PQ{r)/k{r) varies in the same way as Pq{z) in the Cartesian case 
without gravitation. The effects of the different geometry show up as a 
bending of the out streamers towards the equatorial plane. 
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a) 




0,5 1 1,5 2 2,5 3 3,5 4 4,5 5 5,5 

rsin(9) 



b) 

2r 




0,5 1 1,5 2 2,5 3 3,5 4 4,5 5 5,5 

rsin(9) 

Figure 6. Magnetic field fines (contourplot of ^{r,9)) for triple helmet streamers in 
spfierical geometry (see text) 
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In Figure lb) (see Table H case b) we present a solution with cusp 
structure in spherical geometry with the total pressure 



Po{r) 



k{r 



-r+lyn 



+ b 



for r < r, 



k{r)b for r > r^^^, 



cusp 
sp 



(41) 



and use the same parameters as in Figure ^). Again the bending of 
the outer streamers is obvious. 

We remark that prescribing the parameters ci,C2 and pq{Rs) is 
equivalent to prescribing the normal component of the magnetic field at 
the solar surface. This property of the solutions could be used to match 
measured magnetic fields and generate the corresponding streamer 
structure. 



4. Conclusions and outlook 

We have developed a comparatively simple procedure based on the 
method of asymptotic expansion to calculate equilibria of triple helmet 
streamer structures. The procedure represents a generalization of the 
asymptotic expansion method used in the theory of magnetotail equi- 
libria (Schindler, 1972; Birn et al., 1975; Schindler und Birn, 1982). We 
have applied this procedure to the calculation of stationary states of 
triple helmet streamers under various conditions and investigated the 
dependence of the solution properties on the model parameters, the 
influence of the solar gravitational field, the inclusion of field aligned 
flow and on the basic geometry used. 

As a next step we plan to use the calculated configurations as start 
equilibria for MHD simulations. These simulations will answer the ques- 
tion of stability which we did not address in this paper and provide 
important insights into possible mechanisms of eruptive phenomena in 
the solar atmosphere, especially those connected with a triple structure 
of the solar corona. The simulations should have relevance both for huge 
coronal mass ejections and very small but frequent eruptions which may 
represent the source mechanism of the slow solar wind. 
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Appendix 

Here we calculate the terms and ^ which where not explicitly 
given in equation ([iq). 



Middle Streamer 

With c = — Cl, k = ki and xq = one finds: 

dA I f 1 cix [p^ , 

a~ = ~ \ ^ tanh(W —cix^ 

dA 1 



dki kici 

Outer Streamers 

With c = C2, k = k2 and xq = xq2 (see equation (p3|)) one finds: 

dA (2v^\/^ciy^cosh(a;2) + sinh(a2)po\/2^/arY^2;ciC2 - 2 sinh(a2)(ci + C2] 

2(p;;/'^ciyjcosh(a2)c2 



= ( sinh( Q2 ) \/ — C2 + sinh( 02 ) \ — ci 
dk2 V V \ Po 

'p^yfaici cosh(a2)) j [^/po ^fa~\'k.2 Cl cosh( 02)02) 



Therefore we used as abbreviations: 
Po - ^2 exp(c2^s) 



Po 



Q<2 := 7:\/Po ( - V2x ci C2 + 2arctanh (a/oi) \ — C2 
2 V \ Po 

+ 2 arctanh (-^0]") / 
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